Generation of uncorrelated random scale-free networks 
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Uncorrelated random scale-free networks are useful null models to check the accuracy an the 
analytical solutions of dynamical processes defined on complex networks. We propose and analyze 
a model capable to generate random uncorrelated scale-free networks with no multiple and self- 
connections. The model is based on the classical configuration model, with an additional restriction 
on the maximum possible degree of the vertices. We check numerically that the proposed model 
indeed generates scale-free networks with no two and three vertex correlations, as measured by the 
average degree of the nearest neighbors and the clustering coefficient of the vertices of degree fe, 
respectively. 
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Complex networks constitute a general framework for 
the topological characterization of many natural and 
technological systems whose complexity prevents a more 
detailed microscopic description P, U Q • Within this 
framework, these systems are represented in terms of 
networks or graphs 4], in which vertices stand for the 
units composing the system, while edges among vertices 
represent the interactions or relations between pairs of 
units. The focus is thus shift to the topological charac- 
terization of the representative network, a task which is 
largely more feasible and yields, nevertheless, a notice- 
able amount of information on the structure and prop- 
erties of the original system. The empirical analysis of 
many real complex networks has unveiled the presence 
of several typical properties, widely found in systems be- 
longing to a large variety of realms. One of the most 
relevant is given bythe scale-free nature of the degree dis- 
tribution P(fc) 'l',^,'5'|, defined as the probabiHty that a 
randomly chosen vertex has degree k (i.e. it is connected 
to other k vertices). In mathematical terms, the scale- 
free property translates into a power-law function of the 
form 



P{k) - fc-^, 



(1) 



where 7 is a characteristic degree exponent. The presence 
of a scale-free degree distribution can have an important 
impact on the behavior of dynamical processes taking 
place on top of the network. Indeed, scale-free networks 
with exponent 7 in the range 2 < 7 < 3 show large 
fluctuations in their degrees, evident in the presence of 
a diverging second moment, (fc^), in the infinite network 
size limit N 00. This divergency, in turn, show up in a 
remarkable weakness of the network in front of targeted 
attacks 0,0] or the propagation of infectious agents 

It has been recently realized that, besides their degree 
distribution, real networks are also characterized by the 
presence of degree correlations. This translates in the 
observation that the degrees at the end points of any 
given edge are not usually independent. This kind of 
degree-degree correlations can be theoretically expressed 
in terms of the conditional probability P{k'\k) that a 



vertex of degree k is connected to a vertex of degree k' . 
From a numerical point of view, it is more convenient to 
characterize degree-degree correlations by means of the 
average degree of the nearest neighbors of the vertices of 
degree k, which is formally defined as |110|] 



fc„n(fc) =^fc'P(fc'|fc). 



(2) 



Degree-degree correlations has led to a first classifica- 
tion of complex networks according to this property . 
Thus, when A;„„(fc) is an increasing function of fc, the cor- 
responding network is said to exhibit assortative mixing 
by degree, i.e. highly connected vertices are preferentially 
connected to highly connected vertices and vice-versa, 
while a decreasing fc„„(fc) function is typical of disas- 
sortative mixing, being highly connected vertices more 
probably connected to poorly connected ones. For un- 
correlated networks, the degrees at the end points of any 
edge are completely independent. Therefore, the con- 
ditional probability P{k'\k) can be simply estimated as 
the probability that any edge points to a vertex of degree 
fc', leading to Pnc(fc'|fc) = k'P{k')/{k), independent of k. 
Inserting this equation in Eq. the average nearest 
neighbors degree reads 



(fc) ■ 



(3) 



that is, independent of the degree fc. 

Analogously, from a theoretical point of view, correla- 
tions concerning three vertices can be characterized by 
means of the conditional probability P(fc",fc'|fc) that a 
vertex of degree k is simultaneously connected to two 
vertices of degrees fc' and fc". We can estimate this kind 
of three point correlations by means of the clustering co- 
efficient of the vertices of degree fc, c(fc) defined 
as the probability that two neighbors of a vertex of de- 
gree fc are also neighbors themselves. This function can 
be formally written as 



5(fc) = P{k",k'\k)p 



k',k", 



(4) 



k',k" 



2 



where pk',k" is the probabiUty that vertices k' and k" 
are connected given that they have a common neighbor 
[T^ . For random uncorrelated networks, the three vertex 
conditional probability can be written as Pnc(fc", k'\k) — 
Pnc{k"\k)Pnc{k'\k), for k > 1. The connection prob- 
ability can also be computed in this case as pk'.k" = 
[k' — l){k" — l)/{k)N, where the term —1 comes from 
the fact that one of the connections of each vertex has 
already been used [l3 . ll5LIT^ . From the above relations, 
the clustering coefficient becomes 



{k)^N 



(5) 



As in the previous case, for uncorrelated networks, the 
function c{k) is constant and independent of k. There- 
fore, any non trivial dependence of the functions A;„„(/c) 
and c(fc) on the degree is a signature of the presence of 
two and three point correlations respectively. 

While most real networks show indeed the presence 
of correlations, uncorrelated random networks are never- 
theless equally important from a practical point of view, 
especially as null network models in which to test the 
behavior of dynamical systems whose analytic solution 
is usually available only in the absence of correlations 
0,0,11, 01- Therefore, it becomes an interesting issue 
the possibility to generate random networks which have a 
guaranteed lack of correlations. In the particular case of 
scale-free networks, however, finding such algorithms is 
far more difficult than one would expect a priori. In this 
paper, we observe that classical algorithms, which are 
supposed to generate uncorrelated networks, do, indeed, 
generate correlations when the desired degree distribu- 
tion is scale-free and no more than one edge is allowed 
between any two vertices 0,0. To solve this drawback, 
we present and test an algorithm capable to generate un- 
correlated scale-free networks. 

The classical algorithm to construct random networks 
with any prescribed degree distribution P(k) is the so- 
called configuration model (CM) il5|,|2d,|2ll|2l|2l. To 
construct a network with the original definition of this al- 
gorithm, we start assigning to each vertex z, in a set of N 
vertices, a random degree ki drawn from the probability 
distribution P{k), with m < ki < N (no vertex can have 
a degree larger than N) , and imposing the constraint that 
the sum ki must be even. The network is completed 
by randomly connecting the vertices with ki/2 edges, 
respecting the preassigned degrees. The result of this 
construction is a random network whose degrees are, by 
definition, distributed according to P{k) and in which, 
in principle, there are no degree correlations, given the 
random nature of the edge assignment. 

While this prescription works well for bounded degree 
distributions, in which (A;^) is finite, one has to be more 
careful when dealing with networks with a scale-free dis- 
tribution, which, for 2 < 7 < 3, yield diverging fluc- 
tuations, {k'^) — > cx), in the infinite network size limit. 
In fact, it is easy to see that, if the second moment of 
the degree distribution diverges, a completely random 
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FIG. 1: Average nearest neighbor degree of the vertices of 
degree k, k„„{k) (a) and average clustering coefficient c{k) (b) 
for the original CM algorithm with different degree exponents 
7. Network size is TV = IC. 



assignment of edges leads to the construction of an un- 
correlated network, but in which a non-negligible frac- 
tion of self-connections (a vertex joined to itself) and 
multiple connections (two vertices connected by more 
than one edge) are present 0|. While multiple and 
self-connections are completely natural in mathemati- 
cal graph theory j^l, they are somewhat undesired for 
simulation purposes, since most real network do not dis- 
play such structures, and also in order to avoid ambigu- 
ities in the definition of the network and any dynamics 
on top of it. This situation can be avoided by impos- 
ing the additional constraint of forbidding multiple and 
self-connections. This constraint, however, has the nega- 
tive side effect of introducing correlations in the network 
IT^ . As an example of this fact, in Fig. ^ we show 
the functions fc„„(A:) and c(fc) computed from numerical 
simulations of the CM algorithm with no multiple and 
self-connections for different 7 exponents and fixed net- 
work size N = 10^. As we can observe, for 7 > 3, which 
corresponds to an effectively bounded degree distribution 
with finite (fc^), both functions are almost flat, signaling 
an evident lack of correlations. On the other hand, for 
values 7 < 3 there is a clear presence of correlations. This 
correlations have a mixed disassortative nature: vertices 
with many connections tend to be connected to vertices 
with few connections, while low degree vertices connect 
equally with vertices of any degree. 
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FIG. 2: Average nearest neighbor degree of the vertices of 
degree k, k„n{k) (a) and average clustering coefficient c{k) 
(b) for the UCM algorithm with different degree exponents 7. 
Network size is Af = lO''. 



The origin of this phenomenon can be traced back to 
the effects of the cut-off (or maximum expected degree) 
kc{N) in a network of size N. In fact, it is possible to 
show that in order to have no correlations in the absence 
of multiple and self-connections, a scale-free network with 
degree distribution P{k) ~ k^^ and size N must have a 
cut-off scaling at most as ks{N) ~ iV^/^ (the so-called 
structural cut-off) 24]. For a power-law network gener- 
ated using the CM algorithm defined above (i.e. gener- 
ating random degrees in the range m < ki < TV), simple 
extreme value theory arguments show in fact that 



kc{N) ~ N^/i-i-^). 



(6) 



For 7 < 3, we have that kc{N) > N^^^ and therefore it 
is impossible to avoid the presence of correlations. Only 
for the particular case 7 > 3 we recover kc{N) < N^^^, 
which explains the lack of correlations observed in Fig. ^ 
for 7 = 3.5. 

Since it is the maximum possible value of the degrees 
in the network that rules the presence or absence of cor- 
relations in a random network with no multiple or self- 
connections, we propose the following uncorrelated con- 
figuration model (UCM) in order to generate random un- 
correlated scale-free networks: 

I. Assign to each vertex z in a set of A'' initially dis- 
connected vertices a degree ki, extracted from the 



FIG. 3: Numerical average clustering coefficient (c) as a func- 
tion of the corresponding theoretical value, given by Eq. 
for the CM (hollow symbols) and the UCM (full symbols) al- 
gorithms. The different points for each value of 7 represent 
different network sizes iV = 10^, 3 x 10^, lO", 3 x 10*, and 
10^ 



probability distribution P{k) ^ k , and subject 
to the constraints m < ki < N^^^ and ki even. 

2. Construct the network by randomly connecting the 
vertices with '^iki/2 edges, respecting the pre- 
assigned degrees and avoiding multiple and self- 
connections. 

The constraint on the maximum possible degree of the 
vertices ensures that kc{N) ~ N^^^, allowing for the pos- 
sibility to construct uncorrelated networks. As an addi- 
tional numerical optimization of this algorithm, we also 
impose the minimum degree m = 2 to generate connected 
networks with probability one |2^ |2^ . 

In Fig.|2]we check for the presence of correlations in the 
UCM model for scale-free networks. As we can observe, 
both correlation functions show an almost flat behavior 
for all values of the degree exponent 7, compatible with 
the lack of correlations at the two and three vertex levels. 

We have additionally explored the validity of the ex- 
pression for the average clustering coefficient [2^, (c), 
defined as 



(c)=^P(A:)c(fc), 



(7) 



which, for random uncorrelated networks takes the form 
given by Eq. lO 15]. For scale-free networks with a gen- 
eral cut-off kc{N), we have that, in the large A^ limit, 
{k'^) ~ kc{N)^~'' . Therefore, for random networks gen- 
erated with the classical CM model, in which kc(N) ~ 
^1/(7-1)^ we have that (c)cm A^(7-37)/(7-i). This ex- 
pression is clearly anomalous for 7 < 7/3, since it leads 
to a diverging clustering coefficient for large A", while, 
by definition, this magnitude, being a probability, must 
be smaller that one. This anomaly vanishes in the UCM 
prescription. In this case, we have that kc{N) ^ A^/^ 
for any value of 7, leading to (c)ucm ~ N'^~^, which is a 
decreasing function of the network size for any 7 > 2. 



4 



In Fig. 131 we plot the average clustering coefRcient ob- 
tained from numerical simulations of the CM and UCM 
algorithms as a function of the theoretical value, Eq. Q , 
for different values of 7 and different network sizes N. We 
can observe that, while the results for the uncorrelated 
UCM model nicely collapse into the diagonal line in the 
plot, meaning that the numerical values are almost equal 
to their theoretical counterparts, noticeable departures 
are observed for the implicitly correlated CM algorithm. 

To sum up, in this paper we have presented a model to 
generate uncorrelated random networks with no multiple 
and self-connections and arbitrary degree distribution. 
The lack of correlations is especially relevant for the case 
of scale-free networks. In this case, our algorithm is ca- 
pable to generate networks with flat correlation functions 
and an average clustering coefficient in good agreement 
with theoretical predictions. Our algorithm is potentially 
interesting in order to check the accuracy of the many 



analytical solutions of dynamical processes taking place 
on top of complex networks, which are usually found in 
the uncorrelated limit and, which, up to now, lacked a 
proper benchmark to check the results for degree expo- 
nents 7 < 3. 
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